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Abstract 


This  study  uses  linearized  equations  of  motion  for  a  rigid  body  with  an  attached 
spring-mass-damper  to  maximize  the  decay  rate  of  a  satellite’s  coning  motion.  An  analysis 
of  the  numerical  eigenvalues  is  presented  which  leads  to  an  optimal  relationship  between 
relevant  parameters-damper  placement,  spring  constant,  damping  coefficient,  system 
moments  of  inertia,  and  damper  mass  fraction.  The  coupled  system’s  eigenvalues  do  not 
provide  ti'uly  critical  damping,  thus  the  real  eigenvalue  parts  are  minimized  in  order  to 
achieve  damping  which  requires  the  minimum  amount  of  time.  A  comparison  between  this 
optimal  design  method  and  a  cla.ssical  method  concludes  a  noticeable  improvement  in 
damping  perfonuance  for  optimized  systems. 


I.  Introduction 


Background 

Since  the  launch  of  the  first  satellite  in  1957  by  the  Soviet  Union,  dynamicists  have 
investigated  the  motion  of  rigid  bodies  as  they  apply  to  spin-stabilized  satellites.  Most 
satellites  rely  upon  their  ability  to  fix  an  axis  toward  earth  in  order  to  accomplish  such 
missions  as  communications  and  reconnaissance.  A  spin-stabilized  satellite  accomplishes 
this  through  a  large  inherent  angular  momentum,  leaving  perturbing  moments  relatively 
small  and  ineffective.  The  need  for  effective  satellite  attitude  control  was  fuUy  realized 
with  the  failure  of  the  United  States'  Explorer  I. 

The  first  major  tenet  in  the  analysis  of  the  dynamics  of  rotating  bodies  concludes 
that  a  rigid  body  spun  about  an  intermediate  moment  of  inertia  principal  axis  is  unstable 
(17:484).  Adding  energy  dissipation  ensures  that  only  rotation  about  the  principal  axis  of 
maximum  inertia  is  stable  (17:485).  This  was  the  case  of  the  Explorer  I  satellite,  intended 
to  spin  about  its  longitudinal  (minor)  axis.  Long,  flexible  antennas  acted  as  the  energy 
dissipation  mechanism,  causing  a  coning  angle  of  approximately  60  degrees  within  its  first 
orbit  (17:484).  With  this  failure,  designers  turned  their  attention  to  the  design  of  nutation 
dampers  as  a  passive  means  of  improving  attitude  control  about  the  axis  of  maximum 
moment  of  inertia. 

In  modem  practice,  many  satellites  use  advanced  spin-stabilization  methods  such 
as  dual-spin,  momentum  wheels,  or  thrust  control.  Shifts  toward  simpler,  less  costly 
satellites  may  again  utilize  passively  damped  spin-stabilized  systems.  Nutation  dampers. 


similai’  in  design  but  different  in  orientation  to  precession  dampers,  are  devices  which 
increase  energy  dissipation,  tuned  to  return  the  satellite  to  a  state  of  pure  spin  (25:132). 
These  passive  damping  elements  often  resemble  a  spring-mass-damper  system  in  a  fluid- 
filled  tube  and  are  desirable  since  they  do  not  rely  on  complicated  control  systems  or  use 
precious  fuel  (25:133).  Other  damper  types  include  the  viscous  ring  and  the  two  degree  of 
freedom  spherical  pendulum.  The  addition  of  a  nutation  damper  assists  a  spin-stabilized 
satellite  by  damping  out  coning  motions  caused  by  booster  separation  (14:543),  docking, 
crew  motions  (23:1221),  eddy-currents,  or  magnetic  torques  (26:2169).  By  selecting 
proper  components  of  the  nutation  damper,  it  can  efficiently  control  the  attitude  of  a  spin- 
stabilized  satellite.  A  classical  design  method  involves  tuning  the  damper  only  to  the  rigid 
body’s  natural  frequency,  disregarding  coupled  effects  (6:166,  20:298,  21:219,  11:52).  In 
an  uncoupled  model,  the  spring-mass-damper  has  no  effect  upon  the  motion  of  the  rigid 
body,  an  assumption  that  should  be  avoided. 

A  primaiy  area  of  interest  in  the  design  of  nutation  dampers  is  the  relationship 
between  the  spring-mass-damper  parameters.  It  is  the  critical  damping  coefficient  which 
ensures  the  maximum  energy  dissipation  in  the  minimal  amount  of  time.  An  underdaraped 
system  will  eventually  damp  out,  but  will  require  oscillation  of  the  system  over  a  relatively 
long  period  of  time.  Overdamped  systems  will  directly  damp  out,  but  at  a  very  slow  rate. 
Thus,  for  satellite  applications,  it  is  desirable  to  select  a  spring-mass-damper  which  will 
minimize  the  amount  of  time  necessary  to  damp  out  oscillations. 
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History 


The  failure  of  the  Explorer  I  satellite  in  the  late  1950s  may  have  prompted  much  of 
the  work  of  the  following  decade.  The  effectiveness  of  nutation  dampers  in  reducing 
coning  motions  was  well  known,  however  damper  design  was  often  deficient.  In  1962, 
dynamicist  W.  R.  Haseltine  concluded  that  nutation  damper  design  was  “veiy  much  a  cut- 
and-tiy  proposition”  (14:543).  Stability  analyses  led  the  effort,  often  limited  in  scope  by 
computer  capability.  By  1964,  dynamically  stable  systems  were  being  designed,  but  the 
determination  of  the  degree  of  stabilization  was  not  as  straightforward  (24:588).  Due  to 
the  nonlinearities  of  this  problem,  typical  stability  analyses  were  not  exact. 
Approximations,  equation  of  motion  simplifications,  and  computer  simulation  helped 
analysts  study  the  satellite  model. 

Efforts  were  primarily  aimed  at  developing  a  device  which  maximized  the  amount 
of  energy  dissipated  per  pound  of  weight  (21:219).  Several  researchers  studied  the  effects 
of  multiple  dampers  (13:545,  19:848)  and  other  damper  types  (1:456,  2:50,  3:1,  4:383, 
23:1221)  to  include  both  viscous  ring  and  pendulum  dampers.  Viscous  ring  dampers,  in 
the  case  of  the  Telstar  satellite,  incorporate  two  curved  aluminum  tubes  filled  with  neon 
gas  at  one  atmosphere  of  pressure.  Each  tube  contains  a  single  tungsten  ball  which 
dissipates  energy  as  heat  while  moving  through  the  viscous  gas  (26:2182).  The  pendulum 
damper  acts  as  a  combined  precession  and  nutation  damper  due  to  its  two-direction 
degree  of  freedom  (2:50).  In  the  Orbiting  Solar  Observatory,  a  pendulum  damper  is  used 
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which  has  a  circular  steel  arm,  spherical  zinc  bob,  and  silicone  damping  fluid  (11:51).  A 


third  type  of  damper,  the  classic  ball-in-tube,  simplifies  the  analysis  and  has  also  been 
studied  (6:165). 


Problem  Statement 

The  equations  of  motion  of  a  rotating  body  with  an  attached  spring-mass-damper 
have  been  previously  developed  (15:65),  but  much  of  the  analysis  was  limited  to  steady 
state  motion  of  the  system  or  approximate  solutions.  Other  works  have  investigated  the 
ttansient  motion  of  the  system  and  defined  stable  motion  as  functions  of  the  spring's 
stiffness  (7:3,  8:25).  Several  studies  have  proven  that  tuning  the  spring-mass-damper  to 
the  satellite’s  natural  frequency  improves  the  .system  respon.se  (6:165,  11:52,  20:298, 
21:219),  but  few  have  determined  an  optimal  solution. 

This  study  first  reviews  the  derivation  of  the  equations  of  motion  of  a  rigid  body 
with  an  attached  ball-in-tube  damper  in  general  form.  Next,  a  summary  of  the  non- 
dimensionalizations  (8:12)  leaves  the  equations  of  motion  in  a  more  usable  state. 
Lineaiization  of  the  equations  of  motion  about  an  equilibrium  reduces  the  equations  of 
motion  into  their  final  fonn.  Stability  of  the  system  and  its  importance  are  discussed. 

Using  the  linearized  equations  of  motion  of  the  system,  an  analysis  is  conducted  of 
the  eigenvalues  and  eigenvectors.  Two  pairs  of  eigenvalues  for  this  fourth  order  system 
are  used  to  detennine  the  type  and  degree  of  damping  each  parameter  affects.  An 
eigenvector  analysis  verifies  the  high  degree  of  coupling  between  the  spring-mass-daraper 
and  rigid  body,  indicating  that  uncoupled  tuning  methods  are  not  entirely  sufficient. 

After  the  initial  discussion  of  parameters  and  damping,  an  optimization  scheme  is 
introduced  which  endeavors  to  find  the  parameter  set  which  induces  minimum  damping 
time.  The  system’s  real  eigenvalue  parts  are  used  as  the  minimized  parameter  in  order  to 
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induce  the  greatest  system  stability.  Lastly,  this  optimal,  exact  solution  method  is 
compared  to  historical  tuning  techniques  which  assume  a  physically  uncoupled  system. 
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//.  Formulation 


Several  dynamicists  have  developed  equations  of  motion  for  rigid  bodies  with 
attached  spring-mass-dampers.  However,  many  used  approximations  or  made 
assumptions  in  order  to  avoid  nonlinearities  during  the  development  of  the  equations  of 
motion.  This  study  uses  exact  equations  of  motion  for  a  rigid  body  with  an  attached  ball- 
in-tube  spring-mass-damper.  Next,  the  parameters  are  non-dimensionalized  (8:12-15). 
The  resulting  equations  are  highly  nonlinear  and  require  linearization  about  a  desirable 
equilibrium  solution.  The  stability  criteria  are  discussed  and  the  importance  of  designing 
an  inherently  stable  system  is  emphasized.  This  chapter  provides  an  oveiwiew  of  the 
development  of  the  equations  of  motion  used  in  further  analysis. 

Equations  of  Motion 

In  order  to  complete  an  analysis  of  the  system,  the  equations  of  motion  must  first 
be  developed.  An  overview  of  these  relations  is  presented  which  is  most  applicable  to  the 
current  analysis.  The  majority  of  this  derivation  stems  from  Hughes'  text  (15:62-64), 
while  additional  refinement  was  accomplished  by  Chinnery  and  Hall  in  their  paper  (7:2). 
Both  are  recommended  to  the  interested  reader. 

A  schematic  of  the  system  model  is  provided  in  Figure  1,  showing  the  rigid  body, 
^6,  and  attached  point  mass,  (P.  Within  the  spring-mass-damper  (nutation  damper),  the 
point  mass,  nip,  is  allowed  to  move  freely  along  the  axial  direction,  n.  This  motion  is 
governed  by  the  damping  coefficient,  c,  and  the  spring  constant,  k.  The  particle’s  motion 
within  the  spring-mass-damper  provides  the  energy  dissipation  mechanism  necessary  to 
stabilize  the  nominal  spin  about  67. 
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Hughes  develops  the  vectorial  equations  of  motion  by  first  defining  the  linear 

momenta  of  the  rigid  body,  and  the  spring-mass-damper,  pp 
P«  =  j (v  +  CO  X  r)din  =  H-  O)  x 

«  (1) 

P(P  =  (v  -f  0)  X  -1-  in) 

The  total  linear  momentum  of  the  system  is  then 

P  =  P$  +  P(p  -  -  c  X  0)  -t-  m^xn  (2) 

Additionally,  we  are  interested  in  the  component  of  /;^  along  the  damper’s  axial  direction, 
denoted  by 

p„  =  n  ■  pj,  =  Wp  (n  •  V  -  n  X  b  •  00  +  i)  (3) 

The  angular  momentum  is  next  defined  as 

hj  =  jrx(v+a)x r):/m  =  x v -f  J ^  •  o) 

(4) 

h  =  hj+rj,  xp^  =  cxv-i-J(04'  in^xh  x  n 

where  c  and  J  are  the  first  and  second  moments  of  inertia  of  the  system  about  the  origin. 
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The  dimensional  equations  of  motion  as  derived  by  Hughes  (15:64)  are 


p  =  mv„  -  c^co  +  m^xn 
h  =  c^v„  +  Jo)  +  Wpib’^n 

P,i  -  ~ 

(6) 

p  =  -0)  p 
h  =  -o)^h-v/p 

p„  =  (v,,  -  r^^O))  -  cx  -  kx 

Note  that  the  superscript  “x”  denotes  the  skew  symmetric  matrix  of  the  listed  vector.  It  is 
a  simplified  method  for  computing  the  cross  product  through  matrix  multiplication.  For 
example,  the  angular  velocity,  O),  crossed  with  another  vector  can  be  computed  using  O)^ 


CO 


X  _ 


0 

(03 


-0)2 


-CO  3  (0  2 

0  -CO  I 

0),  0 


(7) 


Non-Dimensionalization 

The  equations  of  motion  are  non-dimensionahzed  at  this  point  to  simplify  later 
analysis  (7:2).  The  first  step  is  to  assume  that  there  are  no  external  forces  or  torques,  thus 
linear  and  angular  momenta  are  constant.  The  dimensional  variables  (denoted  by  a 
superscript  star)  as  functions  of  their  dimensionless  counterparts  are  found  on  the 
following  page. 


V 


V 


p  — 

m*p 

♦  ,♦ 

c  /j 

c* 

* 

m 

c  — 

’K  ,  * 

in  h 

b  = 

mb' 

h  _ 

k'll' 

l2 

A.  — 

♦  1  *7 

m  h 

e'  =  1-e 


(9) 

(10) 


The  parameter  b  characterizes  the  distance  from  the  origin  to  the  position  of  the  nutation 
damper  (Figure  1)  and  is  closely  related  to  the  radius  of  gyration,  kg,  about  the  e2  axis 
since  (7:2) 

/2=mX'  (11) 


The  dimensionless  form  of  the  moment  of  inertia  matrix  is 


1 1  +e/«"  -Ebx  0 

J  =  -Ebx  1  0 

0  0  /j  +ebx~ 

and  the  final  non-dimensional  equations  of  motion  are 

1  X 

p  =  v„  — c  0)  +£xn 

h  =  c’^v^  -t- J®  +d?xb^n 
p„  =  e(n^v„  -  -i-i] 

p  =  -®^p 
h  =  -®’^h-v/p 

=  e®  (''o  “  -cx-kx 


(12) 


(13) 
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Without  external  forces,  linear  momentum  is  consei’ved  and  p  =  0  is  assumed 
without  loss  of  generality.  Thus  Vo  and  O)  are  eliminated  from  the  above  equations  and  a 
fifth  order,  non-dimensional  system  is  obtained 


where 


h=:h^K-^L 

p,i  =  ~zlJ  K~^  n'^^'xn  +hY  K~^  L  +  ei'H  j  -  cx  -  kx 

l^  +  n'^h^K-^h 

x  =  S _ 

t'  +  zbn'^h^K~U)^n 


L  =  h  -eMb^n 


K-'  = 


1  zhx 


A 

^bx 


A 

A 


D,  l\ 
0  0 


0 

0 

_1_ 

A 


Dj  =  /,  +  th{z'  -  £.h)x' 
Di  =  /]  +  ££  bx  ~ 

=  /j  +  ££  bx^ 


(14) 


(15) 


Linearization  of  Equations  of  Motion 

To  analyze  the  system  stability,  a  linear  set  of  equations  must  be  available.  Since 
the  equations  of  motion  are  highly  nonlinear,  they  are  linearized  about  a  known, 
equilibrium  solution  corresponding  to  pure  spin  about  the  Q2  axis: 


'flf 

'O' 

Ih 

1 

fb 

= 

0 

Pn 

0 

X 

EQ 

0 

(16) 
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This  solution  provides  a  desirable,  consistent  orientation  of  the  satellite  without  coning; 
other  solutions  may  exist  (7:3).  In  this  attitude,  the  satellite’s  spin  axis  is  aligned  with  the 
angular  momentum  vector. 

As  expected,  symbolic  linearization  using  Mathematica  (28)  led  to  a  5  x  5  Jacobian 
matrix  having  a  zero  row  and  column  for  h2,  representative  of  a  zero  eigenvalue.  The  row 
and  column  were  removed,  resulting  in  the  following  fourth  order  linear  system 


where 


0 

7,-1 


A 

0 


B  ^  B 
0  0 

c  c/j 
~B 

1  _Zl 

B  eB 


0 

-bz 


A 

0 


B  =  b£-e'l^ 


(17) 


(18) 


Equations  17  and  18  are  the  dynamical  relations  u.sed  in  the  following  analyses. 


Stability  Criteria 

The  general  stability  criterion  states  that,  for  a  quasi-rigid  body  with  very  slow 
energy  dissipation,  only  spin  about  the  major  axis  is  asymptotically  stable.  Additional 
stability  criteria  must  also  be  included.  For  this  case,  in  which  the  nutation  damper 
provides  the  energy  dissipation,  the  analysis  becomes  more  complicated.  Linearization 
was  used  to  find  stability  conditions  by  .Saiychev  and  Sazonov  (20:295)  while  Hughes 
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used  linearization  and  the  Routh-Hurwitz  criteria  (15:151).  Chinnery  verified  these 
conditions  using  a  Liapunov  function  (8:24).  Chinnery  and  Hall  defined  the  stability 
criteria  (7:3)  in  tenns  of  the  dimensionless  parameters  (Equations  9,  10)  as 


L. 


I 


(19) 


However,  we  know  from  the  non-dimensional  that  h  =  U  therefore,  the  minimum  stable 
spring  constant  (7:3)  is  found 


k 


cril 


e-b 


(20) 


With  the  requirements  listed  in  Equations  19  and  20,  the  satellite’s  stability  is 
maintained  but  the  degree  of  stabilization  is  the  focus  of  fois  analysis.  The  degree  of 
stability  can  be  defined  by  the  location  of  the  eigenvalues  in  the  complex  plane  (a  more 
stable  system  has  eigenvalues  farther  into  the  left-half  plane).  Correspondingly,  the 
maximum  degree  of  stability  coincides  with  the  minimum  transient  duration  of  the  coning 
motion  (20:295).  Thus,  further  analysis  attempts  to  find  the  set  of  parameters  which 
maximizes  the  stability  of  the  system. 
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IIL  Damping  Effects 


The  amount  of  energy  dissipation  caused  by  a  nutation  damper  can  easily  be 
altered  by  varying  its  characteristics,  namely  the  spring  constant,  damping  coefficient,  and 
particle  mass.  In  this  problem,  however,  the  nutation  damper  is  attached  to  a  rigid  body, 
the  satellite.  This  introduces  additional  parameters  which  affect  the  overall,  coupled 
motion  of  the  system.  The  following  chapter  discusses  the  importance  of  such 
characteristics  as  the  damper  placement,  system  moments  of  inertia,  and  the  ratio  of 
particle  mass  to  system  ma.ss  are  discussed.  It  also  reviews  a  method  of  eigenvalue 
analysis  used  to  study  damping  and  relates  the  method  to  the  peculiarities  of  this  problem. 
Lastly,  this  chapter  analyzes  how  each  parameter  affects  damping  and  the  system’s  overall 
motion. 

Eigenvalue  and  Eigenvector  Analysis 

Given  four  linearized,  homogenous  equations  defined  in  the  previous  chapter 
(Equation  17),  we  must  now  attempt  to  find  the  system’s  eigenvalues  to  determine 
attainable  degrees  of  stability.  Using  the  eigenvalues  to  detennine  the  damping 
characteristics  is  common  in  dynamical  analysis.  This  method,  in  a  mathematically 
uncoupled  system,  provides  two  second  order  problems  and  a  simple  method  for  analyzing 
the  damping  of  separate  components:  the  rigid  body  and  the  spring-mass-damper.  The 
fourth-order  system  presented  hei’e  provides  distinct,  coupled  eigenvalues.  However, 
second  order  analysis  principals  apply. 
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In  a  second  order  problem,  ihe  characteristic  roots  describe  the  damping  of  the 
simple  system.  The  damping  time  index,  as  proposed  by  Borelli  and  Leliakov  (5:345),  is 
defined  as  the  "largest  real  part  of  the  characteristic  roots  of  the  system  model,  lineaiized 
about  the  equilibrium  attitude."  It  is  simply  the  inverse  of  the  non-dimensional  time  to 
damp  out  the  coning  motion  of  the  system.  For  a  damped,  stable  system,  all  eigenvalues 
lie  in  the  left-half  plane  and  minimum  time  damping  occurs  when  the  real  eigenvalue  part  is 
a  minimum. 

Eigenvalues  may  consist  of  a  real  part,  which  describes  the  inverse  damping  time, 
and  a  complex  part,  equal  to  the  damped  frequency.  In  underdamped  second  order 
systems,  a  complex  conjugate  eigenvalue  pair  appears  in  the  complex  plane,  as  seen  in 
Figure  2.  Underdamped  systems  will  eventually  damp  out,  but  at  a  much  slower  rate  than 
desired.  Overdamping  occurs  when  real  and  distinct  eigenvalues  are  present,  leading  to  a 
lengthy  damping  time.  Lastly,  critical  damping  is  the  optimal  solution  to  a  second  order 
damped  system,  where  two  identical,  real  roots  exist.  This  damping  type  will  induce  a 
zero  amplitude  relatively  quickly  as  seen  below. 
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Figure  2,  Damping  Over  Time  and  Eigenvalues  for  Damping  Types 
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Since  the  four-component  Y  vector  is  composed  of  {hi,  hs,  p„,  x),  the  eigenvector 
can  be  used  to  analyze  the  degree  to  which  each  component  plays  a  role  in  the 
coixesponding  eigenvalue  and,  therefore,  the  system  stability.  The  angular  momenta,  hi 
and  hi,  are  somewhat  representative  of  the  rigid  body’s  motion  (denoted  RB,  above), 
while  the  particle’s  position,  x,  and  momentum,  p,,,  are  descriptive  of  the  spring-mass- 
damper’s  motion  (denoted  SMD,  above). 

By  comparing  the  magnitudes  of  the  components  of  the  eigenvectors,  the  degree 
of  coupling  can  be  deteimined.  As  can  be  seen  in  the  eigenvectors  above,  the  magnitude 
of  the  fourth  component,  the  mass  particle  position,  x,  is  dominant.  This  implies  that  the 
system  is  heavily  dependent  upon  x  and  that  the  system  is  highly  coupled  between  the 
spring-mass-damper  and  rigid  body.  For  clarification  purposes,  the  eigenvalues  of  the 
system’s  angular  momentum  will  be  referred  to  henceforth  as  the  rigid  body  pair  while 
those  of  the  particle  mass  and  momentum  will  be  referred  to  as  the  spring-mass-damper 
pair. 


Independent  Parameter  Effects 

The  first  step  in  analyzing  the  effects  of  the  six  parameters  relating  to  this  problem 
is  to  study  their  independent  effects.  Table  1  reviews  the  parameters  which  make  up  the 
linearized  equations  of  motion  (Equations  17,  18). 

Table  1,  List  of  System  Parameters 


c  damping  coefficient 

k  spring  constant 

b  distance  to  damper 

h  moment  of  inertia  about  ei 

h  moment  of  inertia  about  03 

£  mass  fraction 
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Recall  that  e'  =  1  -  e.  All  of  the  parameters  above  are  included  in  this  optimization  study. 

Results  generally  show  that  one  pair  of  the  eigenvalues  can  be  effectively  modified 
by  changing  a  system  parameter,  while  the  other  pair  remains  relatively  constant.  Figure  3 
characterizes  the  movement  of  the  eigenvalues  over  a  range  of  parameters.  This  figure 
shows  that  one  pair  of  eigenvalues,  in  this  ca.se  the  one  relating  to  the  spring-mass- 
daraper,  becomes  real,  and  thus  critical  damping  is  possible.  Unfortunately,  the  other  pair 
of  eigenvalues,  representing  the  rigid  body,  varies  little  over  a  range  of  parameters  and  will 
never  allow  truly  critical  damping. 


Iitiaginaiy 

f.:> 

X  =  Spring-Mas.s-Dainper 

O  =  Rigid  Body  X*^ 

Figure  3,  General  Variation  of  Eigenvalues 


Since  the  rigid  body’s  eigenvalues  never  become  completely  real,  critical  damping  of  either 
component  by  tuning  the  imaginary  parts  of  the  eigenvectors  to  zero,  is  not  a  desirable 
solution  to  this  problem.  In  a  simple,  uncoupled  spring-mass-damper,  the  critical  damping 
factor  is  found  by  equating  the  imaginary  eigenvalue  part  (the  frequency)  to  zero.  This 
relation  reduces  to 

(22) 


17 


While  this  value  drives  the  particle  mass  toward  equilibrium  quickly,  it  causes  the 
damper’s  effect  on  the  rigid  body  to  be  negligible.  Because  of  the  system’s  coupling,  it  is 
not  desirable  to  concentrate  analysis  upon  either  individual  component,  but  to  study  the 
rigid  body  and  spring-mass-damper  stability  simultaneously. 

Note,  in  the  figure  above,  that  at  the  point  where  the  rigid  body  eigenvalues  are  a 
minimum,  the  two  pairs’  paths  cross.  Analysis  of  specific  values  shows  that  the 
eigenvalues  themselves  truly  do  approach  equality  at  optimum,  not  just  their  paths.  This 
implies  that  simply  equating  the  rigid  body  and  spring-mass-damper’s  frequencies 
optimizes  the  damping.  This  technique  is  one  used  in  classical  tuning  practices,  however  it 
assumes  a  physically  uncoupled  model  in  which  the  spring-mass-damper  has  no  effect 
upon  the  rigid  body. 

Analysis  of  the  parameter  effects  of  this  system  begins  with  a  discussion  of  the 
changes  in  the  eigenvalues  over  each  parameter,  while  other  parameters  are  held  constant. 
Next,  hypotheses  on  their  coupled  effects  are  made  and  studied.  Figures  4-8  depict  a 
general  set  of  eigenvalues  over  a  range  of  the  given  parameter.  The  plots  below  are  not  to 
scale  and  are  only  meant  to  be  representative  of  the  effect  of  the  parameter  on  the 
eigenvalues. 

From  Figure  4,  below,  it  can  be  .seen  that  the  real  part  of  each  rigid  body 
eigenvalue  is  decrea.sed  in  the  complex  plane  while  maintaining  a  near-constant  complex 
part  (frequency).  It  is  evident  Irom  this  effect  that  increasing  the  distance  from  the  rigid 
body’s  center  to  the  spring-mass-damper,  b,  promotes  a  faster  damping  rate.  An 
improved  response  occurs  since,  for  the  torque-free  assumptions  of  the  second  chapter. 
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the  spring-mass-damper  applies  the  sole  force  which,  when  coupled  with  a  distance,  b, 
creates  the  restoring  torque.  Based  upon  these  findings,  the  damper  placement  is 
considered  an  independent  parameter  which  has  a  noticeable  effect  on  the  damping  of 
coning  motions. 

3  < - b  < - 0 
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Figure  4,  Effects  of  the  Spring-Mass-Daniper  Placement,  h 
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Figure  5,  below,  helps  to  describe  the  effect  of  the  damping  coefficient  on  the 
system’s  eigenvalues.  The  eigenvalues  resulting  from  an  increased  damping  coefficient 
decrease  to  a  minimum.  This  minimum  occurs  for  a  relatively  small  damping  coefficient, 
as  the  spring-mass-damper’s  frequency  approaches  that  of  the  rigid  body.  Again,  this 
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implies  that  tuning  the  rigid  body  and  the  spring-mass-damper’s  frequencies  optimizes  the 


system. 


Note  that  at  some  critical  value  for  c,  the  spring-mass-damper’s  eigenvalues 
approach,  then  become  real.  It  is  at  this  point  that  the  spring-mass-damper,  by  itself, 
becomes  critically  damped,  while  the  rigid  body  does  not. 


As  discussed  in  the  stability  section  of  Chapter  2,  the  system  wUl  be  unstable  for 
veiy  small  values  of  k,  the  spring  constant.  Figure  6  includes  the  region  where  k  <  kcnt  and 
shows  that,  for  veiy  small  k,  the  spring-mass-damper’s  eigenvalues  begin  in  an  unstable 
region  (to  the  right  of  the  dashed  line),  become  overdamped  (two  distinct,  real 


20 


eigenvalues),  critically  damped  (two  identical,  real  eigenvalues),  then  underdamped  (a 
complex  conjugate  pair).  Underdamping  occurs  as  a  result  of  a  constant,  small  c,  and  a 
relatively  large  k  for  a  given  set  of  parameters.  Thus,  it  appears  that  the  damping 
coefficient  may  have  to  increase  somewhat  for  increasing  spring  constant.  The  degree  of 
proportionality  between  c  and  k  for  optimal  solutions  will  be  discussed  in  Chapter  4. 


Figure  7  shows  the  trend  of  the  rigid  body  eigenvalues  to  become  more  negative  as 
the  moments  of  inertia  increase  in  proportion  to  each  other.  This  implies  that  an 
asymmetric  body  is  more  stable  than  an  axisymmetric  one.  However,  as  h  increases,  the 
real  eigenvalue  parts  far  exceed  those  produced  when  7?  increases.  Recall  that  Ij  is  the 
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moment  of  inertia  about  the  direction  of  the  spring-mass-damper  as  shown  in  the  System 
Model  (Figure  1). 
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With  a  relatively  small  moment  of  inertia  about  the  e.-!  axis,  the  effect  of  the  spring-mass- 
damper  is  greater  because  of  a  lesser  moment  to  overcome.  Since  typical  satellites  are 
fairly  symmetric,  further  analyses  use  values  of  //  and  /?  which  are  similar. 

The  results  of  the  effects  of  the  moments  of  inertia  could  be  used  by  satellite 
designers  to  select  satellite  component  placements  which  are  along  the  Ca  direction,  giving 
rise  to  a  smaller  relative  h.  Allowing  asymmetry  affords  a  much  larger  airay  of 
alternatives  for  the  designer,  however  for  large  initial  coning  angles  the  nutation  damper 


may  not  always  be  the  most  efficient  (21:222). 


Figure  8  describes  the  effect  of  the  mass  fraction  upon  the  eigenvalues  of  the 


system.  As  the  mass  fraction  increases,  the  rigid  body  eigenvalue  decreases,  thus 
improving  the  damper’s  effect.  Meanwhile,  the  spring-mass-damper  eigenvalues  converge 
to  a  point,  then  become  complex,  indicating  underdamping.  As  expected,  for  constant 
damper  placement,  damping  coefficient,  spring  constant,  and  moments  of  inertia,  it  is  more 
desirable  to  have  a  large  mass  fraction.  This  result  supports  the  notion  that  a  greater 
particle  mass  will  have  a  greater  effect  on  the  coning  motion,  while  the  stability  criteria  of 
Equation  20  is  maintained. 


Figure  8,  Effects  of  the  Mass  Fraction,  e 
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IV.  Optimization  Problem 


As  discussed  in  Chapter  3,  it  is  apparent  that  the  eigenvalues  representative  of  the 
rigid  body  will  never  become  real  so  no  truly  critical  damping  exists.  Thus,  an 
optimization  problem  arises  in  which  the  lime  until  the  system  damps  out  is  minimized. 
The  optimization  is  accomplished  by  maximizing  the  magnitude  of  the  real  part  of  the  rigid 
body  eigenvalue.  This  chapter  sets  up  the  optimization  problem,  describes  methods  of 
solution,  and  provides  basic  optimization  results. 

Maximization  of  the  Damping  Rate 

The  non-dimensional  damping  rate  is  simply  the  real  part  of  the  eigenvalue  in 
question.  This  is  true  because,  as  in  the  general  .second  order  case,  the  complex  part  of 
the  eigenvalue  equals  the  damped  frequency  while  the  real  part  is  the  damping  rate.  The 
general  expression  is 

F{t)  =  e'^' {A  cos  Cdt  +  B  sin  m)  (23) 

where  the  eigenvalue  is  equivalently 

a  -t-  ©/  (24) 

In  Equation  23,  the  damping  rate  is  clearly  maximized  as  |  a  |  -^  °o  for  an  exponentially 
decaying  function.  This  result  implies  that  the  magnitude  of  the  real  eigenvalue  part  must 
be  maximized  in  order  to  intensify  the  damping  rate  of  coning  motions.  Recall  that,  for 
stable  solutions,  eigenvalues  lie  in  the  left-half  plane,  ensuring  that  a  <  0. 
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Using  the  second  order  damping  characterizations  of  the  previous  chapter  (Figure 
2),  the  solution  to  this  fourth  order  problem  lies  in  a  study  of  its  eigenvalues.  While  the 
two  pairs  of  eigenvalues  are  not  truly  uncoupled,  they  are  characterized  by  the 
components  of  the  state  vector  which  they  represent.  Recall  that  the  rigid  body  eigenvalue 
pair  is  referred  to  as  the  ones  which  repre.sent  the  eigenvectors  for  the  angular  momentum 
components  and  the  spring-mass-damper  eigenvalues  as  the  ones  which  represent  the  mass 
particle  position  and  linear  momentum. 

At  this  point,  all  of  the  parameters  necessary  to  set  up  an  actual  optimization 
problem  have  been  introduced  and  the  model  has  been  shown  to  be  highly  coupled 
(Chapter  3).  The  problem  is  reduced  to  finding  the  minimum  real  part  of  the  rigid  body 
eigenvalues  from  the  Jacobian  matrix  (Equation  17),  given  the  stability  criteria  as 
described  in  Chapter  2  (Equations  19-21).  The  system  parameters  are  b,  c,  k,  I],  I 3,  and  e; 
however  the  analysis  of  Chapter  3  has  shown  that  only  c  and  k  have  noticeably  dependent 
effects  on  an  optimal  solution. 

Minimization  Methods 

Unfortunately,  symbolic  eigenvalues  are  unattainable  using  Mathematica  (28) 
software.  If  the  symbolic  eigenvalues  could  be  found,  equating  the  partial  derivatives  of 
the  real  part  of  the  rigid  body  eigenvalue  with  zero  would  lead  to  an  analytical  minima. 
Results  lead  to  lengthy,  unusable  equations  and  numerical  optimization  is  necessary. 

MatLab  (27)  has  functions  fniin  and  finins,  which  lead  to  efficient,  numerical 
results.  The  function /mm  is  used  to  find  the  minimum  of  a  function  of  one  variable,  using 
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adjustable  tolerances  and  ranges.  The /mm.?  function  is  veiy  similar,  finding  the  minimum 
of  a  multivariate  function  using  unconstrained  nonlinear  optimization.  In  this  function,  the 
Nelder-Meade  simplex  algorithm  uses  the  direct  search  method  to  compute  minima.  The 
Matlab  manual  provides  additional  information  (18:77-80).  As  an  example,  the  finins 
function  could  be  used  to  find  a  set  {c,  k]  yielding  the  minimum  real  eigenvalue  parts  of 
the  rigid  body  pair,  given  the  rigid  body  characteristics  {b,  e,  /;,  /?}• 

Effects  of  Varying  the  Damping  Coefficient 

Based  upon  the  independent  effect  analysis  of  Chapter  3,  the  damping  coefficient 
which  yields  a  minimum  real  eigenvalue  is  relatively  small.  In  addition,  the  magnitudes  of 
the  damping  coefficient  and  spring  constant  are  somewhat  proportional.  This  section 
verifies  these  hypotheses  and  discusses  their  implications  through  illustrations. 

Prior  analysis  using  the  method  of  multiple  scales  has  postulated  that  the  real 
eigenvalue  parts  are  a  concave  fiinction  of  the  damping  coefficient  (6:165).  Using  the 
exact,  linearized  model  relationships  (Equations  17,  18),  numerical  results  support  this 
effect.  As  seen  in  Figure  9  on  the  following  page,  a  very  small  damping  coefficient  results 
in  real  eigenvalue  parts  of  much  greater  magnitude.  To  make  this  plot,  constant  values  of 
b,  /;,  h,  and  e  are  used,  minimizing  the  real  eigenvalue  parts  while  varying  k.  Note  that 
for  a  larger  damping  coefficient,  no  improvements  in  the  damping  rate  are  attained. 
Further  focus  shifts  to  lightly  damped  regions,  where  values  of  c  ai'e  small. 
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Figure  9,  Real  Eigenvalue  Part  Dependence  Upon  the  Damping  Coefficient,  c 


Based  on  this  analysis,  the  system  must  be  lightly  damped  in  order  to  achieve 
improved  damper  response.  Should  designers  impose  “off-the-shelf’  parts  or  limit  the 


damping  coefficient  in  any  way,  sub-optimal  results  are  expected  and  the  spring  constant 
selection  becomes  even  more  important.  A  better  understanding  of  the  relationship 
between  c,  k,  and  their  combined  effects  on  the  eigenvalues  requires  a  study  of  the  effects 
of  varying  the  spring  constant. 
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Effects  of  Varying  the  Spring  Constant 

Recalling  the  minimum  spring  constant  defined  by  Equation  20  (of  the  order  of 
lO"'*)  and  considering  Figure  10,  it  is  apparent  that  only  relatively  small  values  of  k  will 
result  in  a  minimum  damping  time.  Note  that  the  minimum  of  the  real  part  of  the  spring- 
mass-damper's  eigenvalue  occurs  at  a  slightly  larger  value  than  this  critical  k  (for  the 
values  used  here,  kcru  ~  0.00025). 


Figure  10,  Real  Eigenvalue  Farts  Effects  Hased  Upon  the  Spring  Constant,  k 

In  the  figure  above,  the  minimum  real  eigenvalue  approaches  zero  as  the  spring 
constant  exceeds  0.005.  Recall  that,  in  Figure  9,  the  damping  coefficient  was  also 
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relatively  ineffective  as  it  exceeded  the  same  order  of  magnitude.  This  implies  that  the 
two  are  somewhat  proportional,  at  least  in  regions  beyond  the  optimal  solution.  Figure 
11,  confirms  that  the  damping  coefficient  and  spring  constant  are  almost  linearly  related 
for  regions  beyond  optimal,  and  that  a  minimum  damping  coefficient  is  achieved  at  the 
optimal  solution. 


Figure  11,  Damping  Coefficient  and  Spring  Con.stant  Relation  for  Minimum  Reals 


Effects  of  Varying  the  Location  of  the  Spring-Mass-Damper 

The  placement  of  the  spring-mass-damper  in  relation  to  the  center  of  mass  of  the 
rigid  body  has  a  significant  effect  upon  the  real  part  of  the  system’s  eigenvalues,  a  key  to 
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determining  the  time  required  to  damp  out  coning  motions.  While  the  analysis  of  Chapter 
3  held  the  spring  constant  and  damping  coefficient  constant,  the  optimized  effects  of  the 
damper  location  are  now  illustrated. 


The  non-dimensional  distance  from  the  rigid  body's  center  of  mass  to  the 
equilibrium  point  of  the  spring-mass-damper,  b,  plays  an  even  greater  role  when 
minimizing  the  real  eigenvalue  parts  as  functions  of  c  and  k,  as  seen  in  Figure  12. 
Numerical  analysis  shows  that  as  b  increa.ses,  the  eigenvalues’  real  parts  decrease,  thus 
increasing  the  damping  rate  of  the  coning  motion.  Very  small  values  of  the  eigenvalues’ 
real  part  are  attainable  by  extending  the  distance  to  the  spring-mass-damper.  This  occurs 
as  a  result  of  the  torque  caused  by  the  nutation  damper  and  its  proportional  increase  with 
distance. 


Figure  12,  Minimum  Real  Eigenvalues  for  a  Practical  Range  of  b 
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As  seen  in  Figure  13,  the  magnitude  of  the  real  part  of  the  eigenvalues  not  only 
becomes  relatively  large  as  b  is  increased,  but  a  distinct  minimum  occurs  for  specific 
values  of  the  spring  constant.  Numerical  minimization  below  uses  a  given  value  for  b  over 
a  range  of  k,  then  solves  for  the  damping  coefficients  which  yield  the  minimum  real 
eigenvalue  part.  These  results  verify  that  the  minimum  real  eigenvalues  attainable  for 
typical  satellite  configurations  {b  <  1)  lie  within  the  same  spring  constant  region.  As  the 
distance  to  the  damper  increases,  the  coiresponding  optimum  values  of  k  tend  to  increase 
as  well.  The  large  effect  of  extending  the  position  of  the  spring-mass-damper  could  be 
capitalized  upon  by  designers  through  the  placement  of  dampers  at  the  end  of  extendable 
booms  or  antennas. 
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This  chapter  shows  that  an  optimized  solution  is  obtainable  which  maximizes  the 
damping  rate  of  the  system.  Results  show  that  a  relatively  small  c  and  k  are  desirable  and 
can  provide  a  distinct  minimum  to  the  real  eigenvalue  parts  of  the  rigid  body  mode.  It  has 
also  been  verified  that  faster  damping  rales  are  obtainable  as  the  damper’s  location  is 
extended. 
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V.  Comparisons  to  a  Classical  Tuning  Method 


This  chapter  compares  the  optimal  solution  for  the  linearized  equations  of  motion 
to  a  classical  tuning  method.  This  method  assumes  a  completely  uncoupled  system  for  the 
rigid  body  and  spring-mass-damper,  incorporating  a  mathematical  analysis  to  select 
parameters  (6:165,  11:52,  20:298,  21:219).  Analysis  verifies  that  this  practice  is  capable 
of  producing  favorable  damping  results,  but  is  inferior  to  the  optimization  scheme 
presented  here. 


Classical  Nutation  Damper  Design 

Historically,  dynamicists  have  designed  a  nutation  damper  through  a  method  called 
“tuning”.  Nutation  damper  tuning  equates  the  rigid  body  natural  frequency  to  that  of  the 
physically  decoupled  spring-mass-damper.  The  natural  frequency  of  the  rigid  body,  Oi,,  is 
purely  a  function  of  the  system’s  moments  of  inertia  and  is  found  using  Euler’s  equations 
to  be 


CO.,  = 


1(7, -/,)(/, -73) 

7,73 


(25) 


Next,  the  non-dimensional  damped  frequency,  O),/,  is  found  for  the  decoupled  spring-mass- 


damper 


CO, 


(26) 
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After  equating  them,  it  is  a  fairly  simple  calculation  to  detennine  a  set  of  spring  constants 
and  damping  coefficients  which  satisfy  the  tuning  condition  using  the  relation 


c  =  2e^--(0,;-  (27) 

In  Figure  14,  it  appears  that  the  eigenvalues  of  the  rigid  body  and  spring-mass- 
damper  approach  equality  at  the  optimal  solution.  This  finding  supports  tuning  the 
frequencies  of  the  decoupled  rigid  body  and  spring-mass-damper.  The  dashed  line  below 
is  equivalent  to  the  uncoupled  rigid  body  natural  frequency  found  by  Equation  25.  For 
this  analysis,  however,  that  assumption  will  not  be  made  and  it  will  be  shown  that  much 
better  results  are  possible. 
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Actual  Results 


Using  the  physically  decoupled  tuning  method  relationship  for  k  and  c  listed  above, 
damping  responses  of  the  numerical  optimization  and  classical  tuning  are  compared. 
Figure  15  plots  the  relationship  between  the  damping  coefficient  and  spring  constant  for 
optimal  methods,  as  in  Figure  11,  and  classical  tuning.  The  two  methods  result  in  different 
relations  for  c  and  k,  while  both  yield  similar  results  at  the  optimized  solution  (when  c  is  a 
minimum).  Unfortunately,  exact  values  of  the  parameters  are  not  available  at  this  point 


using  the  decoupled  method  alone,  leading  the  designer  to  suboptimal  results. 


Figure  15,  Spring  Constant/Damping  CoeJ'ficient  Relationships 
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A  major  shortcoming  of  the  decoupled  method  lies  in  solutions  when  the  resulting 
damping  coefficient  is  imaginaiy.  This  occurs  when  the  radical  teirn  in  Equation  27  is 
negative,  when 

-<a;-  (28) 

e 

This  region  includes  unstable  values  of  the  spring  constant  but  also  a  noteworthy  region  of 
viable  values  for  k.  For  spring  constant  values  outside  of  this  region,  no  real  value  for  c  is 
obtainable  when  using  the  decoupled  tuning  approach. 

Finally,  a  direct  comparison  between  the  two  tuning  methods  is  made.  Figure  16 
shows  that  classical  tuning  results  in  a  minimum  of  real  eigenvalues,  however  limited  by 
Equation  28.  In  the  plot  below,  a  value  of  b  =  1.0  is  used  to  illustrate  the  trend  in  the 
differences.  The  tuned  curve  uses  values  of  c  which  are  provided  thi'ough  optimization 
results  in  order  to  emphasize  the  improved  perfonnance  of  optimal  methods.  In  practical 
applications,  designers  must  arbitrarily  select  a  value  for  one  parameter  (c  or  k),  then  solve 


for  the  other  using  Equation  27. 
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The  minimum  point  on  the  decoupled  tuning  curve  above  occurs  when  c  is  slightly  greater 
than  zero  (Figures  15  and  16).  This  implies  that  designers  should  select  the  minimum 
damping  coefficient  available.  While  this  may  seem  plausible,  available  hardware  design 
may  limit  its  practicality.  By  using  optimal  tuning,  a  relationship  between  the  damping 
coefficient  and  spring  constant  is  available  which  offers  exact  values  for  c  and  k. 

The  differences  in  the  real  eigenvalue  parts  are  the  key  to  the  realizing  the  benefits 
of  optimal  parameter  selection.  Figure  17  plots  the  net  gain  in  damping  rates  over  the 
same  range  of  the  spring  constant  as  seen  in  Figure  16.  The  net  gain  is  characterized  by 


NetGain  =  ||g  -  fa  (29) 
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Here,  the  function  is  optimized  around  the  damping  coefficient  while  vaiying  the  spring 
constant.  The  primary  benefit  offered  by  optimal  methods  is  an  increase  in  the  optimum 
region.  By  inspection  of  the  classical  tuning  results  (Figures  16  and  17),  a  small  change  in 
k  has  a  shaip  effect  on  the  magnitude  of  the  damping  rate  because  of  the  slope  of  the 
curve  near  the  minima.  Using  optimal  methods,  this  effect  is  reduced  by  providing  a  larger 
region  of  improved  damping  rate  perfonnance. 


In  addition  to  the  benefits  obtained  by  optimal  tuning  methods  for  a  constant 
damper  placement,  even  greater  effects  are  achieved  with  increasing  b.  For  the  same 
range  of  values  for  b  as  discussed  previously.  Figure  1 8  shows  a  clear  distinction  between 
the  damping  rates  obtained  using  optimization.  Since  an  exact  solution  to  c  is  not 
practically  obtainable  using  classical  methods,  a  constant  small  c  is  used  (c  =  0.001).  As 
seen,  this  seriously  limits  the  performance  of  the  damper  when  b  is  less  than  one. 


Figure  18,  Response  of  Real  Eigenvalue  Farts  for  Tuning  Methods 
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Figure  19  provides  a  comparison  of  damping  coefficients  and  spring  constants  for  the 
same  range  of  b.  Note  that  the  uncoupled  method  does  not  allow  c  and  k  to  vary  as  Z?  is 
increased. 


Figure  19,  Damping  Coefficient  and  Spring  Constant  Comparison  for  Tuning  Methods 


These  findings  show  that,  although  the  classical  tuning  method  provides  improved 
system  damping  response  to  coning  motions,  optimization  can  achieve  much  better  results. 
Using  exact,  linear  equations  of  motion,  optimization  schemes  are  created  which  are  useful 
over  all  stable  regions  of  the  parameters  given.  On  the  contrary,  classical  tuning  methods 
are  not  applicable  for  small  values  of  the  spring  constant,  k.  Due  to  an  extreme  sensitivity 
of  the  damping  rate  to  the  spring  constant  in  classical  methods,  any  slight  deviations  in 
selection  of  the  spring  constant  through  Equation  27  could  lead  to  unsatisfactory  results. 
These  slight  discrepancies  come  as  a  result  of  the  limitations  in  selecting  an  arbitrary,  small 
damping  coefficient.  By  using  optimization  schemes,  the  nutation  damper  designer  can 
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obtain  exact  values  for  the  damping  coefficient  and  spring  constant,  limited  only  by  the 
optimization  algorithm  itself. 

Optimization  of  the  exact,  linearized  equations  of  motion  allows  the  satellite 
nutation  damper  designer  to  obtain  values  for  the  system  parameters  which  maximize  the 
rate  of  damping.  The  classical  tuning  method,  in  which  the  rigid  body  and  spring-mass- 
damper  are  decoupled,  provides  a  parameter  set  yielding  inferior  results.  This  chapter 
shows  that  optimization  provides  results  which  are  more  precise  and  promote  faster 
damping  rates. 
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VI.  Conclusion 


Designers  have  realized  the  need  for  good  nutation  damper  designs  over  the  course 
of  spacecraft  evolution.  Efficient  passive  attitude  control  mechanisms  are  critical  to 
providing  simpler,  less  costly  mission  support.  Spin-stabilized  satellites  in  particular  are 
key  prospects  for  the  incoiporation  of  nutation  dampers  to  control  coning  motion.  This 
satellite  “wobble”  may  be  caused  by  external  perturbations  such  as  crew  movement, 
booster  separation,  docking,  or  magnetic  torques.  A  nutation  damper,  acting  as  the 
energy  dissipation  device,  returns  the  satellite’s  angular  momentum  vector  to  an  alignment 
with  its  spin  axis.  In  doing  so,  the  satellite  remains  pointed  to  accomplish  its  designed 
mission. 

The  study  of  a  rigid  body  with  an  attached  spring-mass-damper  has  often  coincided 
with  satellite  research.  While  many  of  the  stability  analyses  of  these  systems  were 
developed  during  the  1960s,  designs  eventually  shifted  toward  more  complicated,  costly 
control  schemes.  Historically,  only  approximate  solutions  have  existed  in  discussions  of 
nutation  dampers  and  studies  have  been  limited  to  estimates  on  their  designs.  In  classical 
damper  design,  the  mathematically  uncoupled  rigid  body  frequency  is  equated,  or  tuned,  to 
that  of  the  spring-mass-damper.  While  this  method  has  been  shown  to  improve  the  cone 
angle  decay  response,  it  is  not  truly  an  optimal  solution. 

With  the  non-dimensional  equations  of  motion,  an  analysis  of  the  eigenvalues  leads 
to  a  stability  analysis  of  the  system.  Based  upon  the  eigenvalue  response,  the  spring-mass- 
damper  can  damp  out  critically  while  the  rigid  body  cannot.  To  achieve  favorable 
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damping  conditions,  the  damping  rate  is  optimized  by  maximizing  the  magnitude  of  the 
real  part  of  the  stable,  rigid  body  eigenvalues. 

A  study  of  individual  parameters  which  describe  the  symmetry,  h  and  h,  the  mass 
fraction,  e,  damper  placement,  b,  damping  coefficient,  c,  and  spring  constant,  k,  illustrates 
their  independent  effects.  The  moments  of  inertia,  h  and  h,  have  a  noticeable  effect  on 
the  performance  of  the  nutation  damper.  A  relatively  asymmetric  rigid  body  provides  a 
greater  damping  rate,  and  when  7;  >  h,  a  much  improved  damping  rate  is  seen.  This 
phenomenon  is  expected  since  the  nutation  damper  studied  has  a  single  degree  of  freedom. 
When  increased,  the  mass  fraction,  e,  provides  a  greater  damping  response  to  coning  as 
well.  As  the  damper  placement,  h,  increases,  the  damping  rate  also  increases  due  to  the 
torque  applied  to  the  rigid  body.  Since  the  effects  of  /;,  /?,  8,  and  b  appear  to  be 
independent,  they  are  fixed  in  the  optimization  function,  leaving  only  the  damping 
coefficient,  c,  and  spring  constant,  k,  as  variable  parameters.  Results  show  that  for  small  c 
and  k,  the  damper  provides  favorable  response  to  coning.  As  the  damping  increases,  the 
system  becomes  overdamped  and  when  the  spring  constant  increases,  the  system  becomes 
underdamped.  This  result  illustrates  that  the  damping  coefficient  and  spring  constant  are 
somewhat  proportional. 

Optimization  reduces  to  maximizing  the  damping  rate  by  minimizing  the  real 
eigenvalue  parts  of  the  rigid  body.  Only  c  and  k  are  allowed  as  optimization  variables, 
although  the  effects  of  b  are  depicted.  It  is  seen  that  optimal  damping  response  occurs 
only  for  small  values  of  c  and  k,  while  verifying  that  the  two  are  almost  lineaiiy  related  for 
non-optimal  regions. 
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By  comparing  the  rigid  body  and  spring-mass-damper  eigenvalues  at  optimum,  the 
two  pairs  approach  equality.  Initially,  classical  tuning  methods,  which  equate  the  two 
frequencies,  seem  to  be  correct.  However,  it  is  clear  that  an  optimized  solution  is  capable 
of  producing  a  remarkable  improvement  in  damping  rates.  In  classical  methods,  a  slight 
deviation  in  parameter  selection  leads  to  marked  differences  in  performance.  On  the 
contrary,  optimal  solutions  offer  an  exact  result  which  yield  the  greatest  attainable 
damping  rate  for  the  decay  of  coning.  The  solutions  pre.sented  here  can  be  modified  to  fit 
other  satellite  models  and  may  provide  satellite  attitude  control  system  designers  with 
improved  perfonnance  of  simple,  passive  nutation  dampers  attached  to  spin-stabilized 
satellites. 
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Appendix  A:  Sample  Mat  Lab  Code 
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%%%%%%%%%%%%%%%%%%%%%Wo%%%%%%Wc%%%%%%%%%%%%%%%%%%%%% 


%  % 

%  This  code  defines  conshints,  tlien  solves  for  values  of  c  which  % 

%  unpart  tlie  minimum  real  eigenvalue  pmt  over  a  nmge  of  k  =  % 

%  {0.002,  0.004}.  It  must  be  used  in  conjunction  witli  the  % 

%  function  F(c).  % 

%  % 

%  The  optunization  tolerance  of  fmin  may  be  defined  as  10  '“  by:  % 

%  fmin(‘F’,emin,cmax,[0,10*(-10)])  % 

%  % 

%  Similai'ly,  tlie  fmins  function  can  be  used  to  solve  for  muUiple  % 

%  variables  by  invoking;  % 

%  fmins(‘MinFn’,[guessl,gue.ss2],[0,]0''(-10)])  % 

%  % 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%e  =  e,  ep  =  e' 

e=0.01; 

ep=l-e; 

13=0.7; 

11=0.6; 
b  =  3; 

B  =  b*e  -  ep^'Ll; 
j  =  0; 

%cmin  and  cmax  define  range  of  search 

cmin  =  0; 

cmax=0.02; 


Vais  =  zeros(100,2); 

for  k  =  0.002:0.00002:0.004 

j=j+i; 

%Find  c  which  yields  minimum  real  eigenvalue  (=  ans) 
fmin('F',cmin,cmax); 

Min=F(ans); 

Vals(j,:)  =  [k  Min]; 
end 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 

%  This  function  defines  ‘F’  as  a  routine  which  will  return  tire 
%  value  for  c  which  results  in  tlie  miniininn  real  eigenvalue  par  t. 

%  It  must  be  used  in  conjunction  with  a  statement  resembling: 

%  fmin(‘F’,cmin,cmax) 

% 


%  Where  [cmin.cmax]  =  range  of  acceptable  values  for  c. 


function  minev  =  F(c) 
global  k  b  e  ep  13 II; 

0; 

-(b*e/Il); 
-(b*e^2/II)-k; 
0]; 


M=  [  0 

(-1+I1)/II 

-(e/Il) 

0 


- 1  +ep/(-(b*e)+ep*I3) 
0 

c/(-(b*e)+ep*I3) 

-(l/(-(b*e)+ep*I3)) 


-(b/(-(b*e)+ep*I3)) 

0 

-c*I3/(e*(-(b*e)+ep*I3)) 

13/(c*(-(b*e)+ep*I3)) 


evals  =  eig(M)'; 

[minev, imin]  =  min(real(evals(:,l))); 
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%  % 

%  This  function  defines  ‘MinFn’  as  a  routine  which  will  return  % 

%  values  for  c  and  k  which  yield  tlie  mlniinuin  real  eigenvtdue  pait.  % 

%  It  may  be  used  in  conjunction  wiili  the  statement  % 

%  fmins(‘MinFn',[guessI,  guess2j)  % 

%  % 

%  Where  [guessl,guess2]  =  initial  estimates  of  c  tuid  k.  % 

%  % 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


function  minev  =  MinFn(v) 
global  e  ep  13  II  b; 
c=v(l); 
k=v(2); 

0; 

-(b*e/Il); 
-(b*e''2/Il)-k; 
0]; 

evals  =  eig(M)'; 

[minev, imin]  =  min(real(eva]s(:,l))); 


M  =  [  0  -l+ep/(-(b*e)+ep*I3)  -(b/(-(b*e)+ep*I3)) 

(-1+I1)/I1  0  0 

-(e/Il)  c/(-(b*e)+ep*I3)  -c*I3/(e*(-(b*e)+ep*I3)) 

0  -(l/(-(b*e)+ep*I3))  I3/(c*(-(b*e)+ep*I3)) 
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Appendix  B:  Supportive  Plots 
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Minimum  Attainable  Real  Eigenvalue  Parts 


Effect  of  Increasing  b  on  the  magnitude  of  the  damping  rate 
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Min  Real  Eigenvalues 


Effect  of  Increasing  c  to  0.008 


Spring  Constant,  k 


50 


Min  Real  Eigenvalues 


5 


Min  Real  Eigenvalues 


Effect  of  Increasing  the  Mass  Fraction 
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